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The Wigner function formalism which relates source size parameters 
to experimental "HBT radii" extracted from two-particle Bose-Einstein 
correlations is generalized to azimuthally deformed and longitudinally tilted 
sources. It is explained how this can be used to complement anisotropic 
flow measurements with relevant space-time information on the source. 
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1. Source size parameters and HBT radii 

Two-particle interferometry, which exploits the Bose-Einstein symmetriza- 
tion effects on the production cross section for pairs of identical bosons, has 
become a powerful tool to extract detailed space-time information about the 
freeze-out configuration of the hot and dense fireballs formed in relativistic 
heavy- ion collisions . Recent exciting data on transverse flow anisotropies 
in non-central heavy-ion collisions and their interpretation as evidence for 
early and efficient thermalization at RHIC energies [2] have generated new 
interest in a better understanding of the space-time structure of the de- 
formed sources created in these collisions. This requires the generalization 
of the HBT interferometry tool from azimuthally symmetric sources (re- 
viewed in [I]) to deformed situations. Applications of this new formalism 
to recent RHIC data are reported in the following article by M. Lisa. 

For chaotic sources (i.e. independent particle emission), the two-particle 
Bose-Einstein correlation function C(p\,p2) = C(q,K) can be expressed 
through a Fourier transform of the emission function S(x, K) which de- 
scribes the single-particle phase-space distribution at freeze-out [TJ. K is 
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the average momentum of the pair while q is the relative momentum be- 
tween the two particles. The Fourier transform in q is restricted by the 
mass-shell constraint q° = E\—Ei = (3 ■ q (where (3 = K/K° is the pair ve- 
locity) and therefore not fully invertible. The extraction of the emission 
function S(x,K) from the measured correlation function C(K,q) thus re- 
quires additional theoretical input pQ. 

If the source is dominated by a single length scale (its "size" ) , the emis- 
sion function can, for every momentum K, be approximated by a Gaussian 
in space-time whose width parameters form a symmetric tensor, the spa- 
tial correlation tensor S^ U {K) = [(x^x u ) — (x^)(x u )](K) = (x^x u ). The 
correlation function is then also a Gaussian in the relative momentum q, 
C(q, K) = \ + exp[— J2i,j=o,s,l HQjRiji-^-)]) where I denotes the beam di- 
rection, o the transverse emission direction K± of the pair, and s the third 
Cartesian direction perpendicular to I and o. The transverse emission direc- 
tion K is characterized by an azimuthal angle with respect to the reaction 
plane formed by the beam axis and the impact parameter b. In collisions 
between spherical nuclei the emission function is reflection symmetric with 
respect to the reaction plane. This and other, more specific symmetries 
of the emission function are most easily expressed in a reaction-plane-fixed 
coordinate system where z points along the beam direction, x along the 
impact parameter, and y perpendicular to the reaction plane. On the other 
hand, the correlation radii Rfj(K) are more easily interpreted in the (o, s, I) 
system because the ^-dependence, being transverse to the pair velocity /3, 
is not affected by the mixing of space and time induced by the mass-shell 
constraint into the Fourier transform. The (o, s, I) and (x, y, z) systems are 
rotated with respect to each other by the azimuthal emission angle <1>. 

The Fourier transform between the emission and correlation functions 
the leads to the following relations between the "HBT radii" Rfj(K) and 
the components of the spatial correlation tensor S^^K) [3]: 

R 2 S = \{S XX + Syy) -\{S XX -Syy) COS{2$) - S X ySlll(2$) 

Rl = \{S xx +S yy ) + \{S xx -S yy )cos{2®) + S xy sni{2$) 

-2f3 x (S tx cos $+S ty sin $) + 0{S tt , 
Rls = s xy cos(2$) - \ (S xx -S yy ) sin(2$) + (3±(S tx sin $-S ty cos $), 
Rf = S zz — 2fiiStz + pfStt, 

R 2 oi = (Sxz-fiiStx) cos $ + (S yz -f3iS ty ) sin® - f3±S tz + (3if3±S tt , 

R 2 sl = (Syt-Pi^coaQ-iSM-ftStJsm*. (1) 

These relations display the explicit ^-dependence arising from the men- 
tioned rotation between the (x, y, z) and (o, s, I) systems, but hide the im- 
plicit ^-dependence of the spatial correlation tensor components S^ V {K) = 
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Sui/(Y, K±, $). The total emission angle dependence of the HBT radii re- 
sults from a combination of both explicit and implicit ^-dependences 013]. 

The implicit ^-dependence of the spatial correlation tensor is restricted 
by symmetries of the source [I]. It is a relativistic effect associated with 
an azimuthal spatial source deformation superimposed by strong transverse 
collective flow which vanishes with the 4 th power of the transverse flow 
velocity vt/c for weak or no collective expansion 

2. Azimuthal oscillations and harmonic analysis 

A full analysis of the symmetry constraints on S flu (Y, K±, $) for symmetric 
collisions between spherical nuclei and for pairs detected in a symmetric 
rapidity window around Y = can be found in Ref . [3] . One finds the 
following most general form for the azimuthal oscillations of the HBT radii: 

R l = R lo + 2 E„=2,4,6,...^,n COs(n$), R? os = 2 T, n =2,4,6,... R ls,n sin(n$), 
R l = R l,o + 2 E„=2,4,6,...^,n cos(ra$), R 2 ol = 2 £n=l,3,5,...-Rrf,n cos(n$), 
Rl = Rlo + 2E„= 2 ,4,6,...^ 2 nCos(ncI>), R 2 sl =2E n=1 ,3 A ...4,r ! ™(^).(2) 

We see that only even or odd sine or cosine terms occur, but no mixtures of 
such terms. Statistical errors in the resolution of the reaction plane angle 
as well as finite angular bin sizes in $ tend to reduce the actually mea- 
sured oscillation amplitudes; fortunately, these dilution effects can be fully 
corrected by a model- independent correction algorithm [2]. A Gaussian 
fit to the thus corrected correlation function, binned in Y, K± and emis- 
sion angle <I>, then yields the "true" HBT radius parameters R?j(Y,K±,$) 
from which the n th order azimuthal oscillation amplitudes are extracted 
via Rfj^K^) = ^E^i2? (F ) ir_ L ,$ i )osc(n$ i ). Here n bin indicates 
the number of (equally spaced) $ bins in the data and osc(n<I>j) stands 
for sin(n<3?j) or cos(n$j) as appropriate, see Eqs. ©. (Note that Nyquist's 
theorem limits the number of harmonics that can be extracted to n < nbin-) 

3. HBT oscillation amplitudes and source shape 

We would like to relate the azimuthal oscillation amplitudes of the 6 HBT 
radius parameters to the geometric and dynamical anisotropies of the source, 
as reflected in the azimuthal oscillations of the 10 independent components 
of the spatial correlation tensor. Their allowed oscillation patterns at midra- 
pidity Y = are given by [1] 

A(*) = \{x 2 +y 2 ) = A +2E„>2,evenA„COs(n$), 
B{<5>) = \(x 2 -y 2 ) = £ +2£n>2,even£nCOs(n$), 
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C(d>) = 


(xy) 


= 2^ n > 2even C n sin(n$), 




(?) 


= L>o+2En>2,evcn- D nCOs(n$) 


= 




2En>l,odd £: nCOs(n$), 


F(*) = 




= 2 En>l,odd i? nS in ( n$ ) 1 


G(*) = 




= 2 En>l,odd G! nCOs(n$), 


#(<!>) = 


(xz) 


= ff +2E n >2,even- ff nCOs(n*), 


/(*) = 


(P) 


= 2E„>2,even / nCOs(n$), 


J(*) = 


(i 2 ) 


= Jo+2En>2,even^nCO S (n$). 



(3) 

The missing terms in the sums over n have amplitudes which are odd func- 
tions of Y and vanish at midrapidity. They do, however, contribute to 
the HBT radii if the data are averaged over a finite, symmetric rapidity 
window around Y = 4 . Their contributions can be eliminated by vary- 
ing the width AY of this rapidity window and extrapolating quadratically 
to Ay — ► 0. Note that Co = Eq = Fq = Gq = Iq = by symmetry, i.e. the 
corresponding components of oscillate around zero. 

The oscillation amplitudes of the HBT radii relate to the oscillation 
amplitudes of the source parameters as follows: For the odd harmonics 
n = 1, 3, 5, . . . we have 

RoLn = ^(Hn-l+H n+ i—I n -i+I n+ i—l3i(E n _i+E n+ i — F n ^i+F n+ i)) 

-(P±G n -PiD n ), (4) 

Rsl,n = ^( — Hn-l+Hn+l+In-l+In+1 — (3l(—E n -i+E n+ i+F n _i+F n+ i)), 

whereas the even harmonics n = 0, 2, 4, . . . satisfy 

Rs,n = (An) + \ ( — B n -2— B n+ 2+C n -2 — C n+ 2) i 
Ro,n = (A n ) + \(B n -2+Bn+2 — C n -2+C n+ 2) 

—(3±(E n -i+E n+ i—F n -i+F n +i) + /3j_{D n ), 

Ros,n = \{~ B n -2+B n+ 2+C n -2+C n+ 2 + (3±(E n _i — E n+ i~ F n _i — F n+ i)), 
R\ n = (J n )-2((3 l G n ) + (pfD n ). (5) 

In these relations it is understood that all negative harmonic coefficients 
n < as well as Co, Eq, Fq, Go and Iq are zero. The angular brackets (. . .) 
indicate an average over a finite, symmetric rapidity window around Y = 0. 
The terms involving the longitudinal pair velocity j3i vanish quadratically 
as the width Ay of that window shrinks to zero. 

Even after extrapolating to Y = in this way, we have still many more 
source parameters than measurable HBT amplitudes. One counts easily 
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that up to n = 2 there are 9 measurable Fourier coefficients which (at Y = 0) 
depend on 19 source amplitudes. From there on, increasing n by 2 yields 
6 additional measured amplitudes which depend on 10 additional source 
amplitudes. This lack of analysis power is an intrinsic weakness of the 
HBT microscope and due to the fundamental restrictions arising from the 
mass-shell constraint q° = (3 ■ q. The reconstruction of the source thus must 
necessarily rely on additional assumptions. 

One such assumption which may not be too unreasonable is that the 
emission duration D = {i 2 ) is approximately independent of emission an- 
gle and that the source is sufficiently smooth that higher order harmonics 
n > 3 of S^y can be neglected. Such source properties would result in the 
"Wiedemann sum rule" [H] 



which can be experimentally tested. If verified for all K± it would provide 
strong support for the underlying assumptions on the source. In this case we 
can measure 3 azimuthally averaged HBT radii and 5 independent oscillation 
amplitudes with n < 2, depending on 14 source parameters of which 5 can 
be eliminated by going to K± = (3± = (see 0] for explicit expressions). 
This makes the geometry of the effective source for particles with K± = 
"almost solvable" |7] , as confirmed by hydrodynamical calculations 9 which 
show that at K± = the effective emission region closely tracks the overall 
geometry of the source even if it is strongly and anisotropically expanding. 

The source geometry can be completely reconstructed from HBT data 
if transverse flow is so weak that all implicit ^-dependence (i.e. all higher 
harmonics n > 1) of can be neglected. In this case one obtains at Y = 
the "geometric relations" [H] 



Aq describes the average transverse size and -Bo (which generates a second- 
order harmic in the transverse HBT radii) the transverse deformation of the 
source. Hq generates a first-order harmonic in the ol and si cross terms and 
describes a longitudinal tilt of the source away from the beam direction jSJ. 
Such a tilt was found in Au+ Au collisions at the AGS ^U] . Its sign yielded 
important information on the kinetic pion production mechanism |Sl llUj. 



(6) 




R ll,l = ~ R sl,l = 2 H = \{%z)o, 

lB = l(x 2 -y 2 } . 



(7) 



4. Conclusions 



Azimuthally sensitive HBT interfrometry is a powerful tool for analyzing 
the dynamic origin and space-time manifestations of the strong anisotropic 
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collective flow seen in single-particle spectra at RHIC. Symmetries strongly 
constrain the azimuthal Fourier series of the emission angle dependent HBT 
radii. This is helpful but not sufficient for a fully model independent recon- 
struction of the source deformations from such data. To separate temporal 
from geometric contributions to the HBT radii, the source must satisfy 
special properties. One can test for them experimentally, albeit not in a 
fully model-independent way. Hydrodynamic simulations 9 (which were 
reported at the meeting but are not included here) show that at RHIC and 
LHC energies emission from the source at non-zero transverse momentum 
is very strongly surface dominated ( "source opacity" ) . To obtain an unam- 
biguous estimate of the transverse deformation of the source at freeze-out 
one should study the ^-oscillations of the HBT radii at very small K±. 
A possible longitudinal tilt of the source away from the beam direction 
manifests itself through first-order harmonics in B?, and R%. The in- 
dependence of the oscillation patterns of the transverse HBT radius param- 
eters can be analyzed to obtain evidence for or against a faster expansion of 
the source in-plane than out-of-plane due to anisotropic collective flow (§]. 
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